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ABSTRACT 



We present properties of the peaks (maxima) of the microwave background 
anisotropies expected in flat and open cold dark matter models. We obtain 
analytical expressions of several topological descriptors: mean number of maxima 
and the probabilty distribution of the gaussian curvature and the eccentricity of 
the peaks. These quantities are calculated as functions of the radiation power 
spectrum, assuming a gaussian distribution of temperature anisotropies. We 
present results for angular resolutions ranging from 5' to 20' (antenna FWHM), 
scales that are relevant for the MAP and COBRAS/SAMBA space missions 
and the ground-based interferometer experiments. Our analysis also includes 
the effects of noise. We find that the number of peaks can discriminate between 
standard CDM models, and that the gaussian curvature distribution provides a 
useful test for these various models, whereas the eccentricity distribution can 
not distinguish between them. 



Subject headings: cosmology: cosmic microwave background — anisotropies: 
peaks 
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1. Introduction 

The standard way to study the microwave background anisotropies (CMB) is based 
on the computation of the radiation power spectrum, i.e., the C^s. The texture of the 
CMB offers a useful alternative to this approach and can be used to test models of galaxy 
formation. Let us consider the excursions of a gaussian random field above a certain 
threshold v — (AT)/(AT)rms- It is of interest to remark that once v is fixed, all the 
topological quantities we will calculate are normalization-independent. Earlier work on the 
properties of peaks in one-dimensional scans and 2D maps of the CMB are due to Zabotin 
& Nasel'skii (1985) and Sazhin (1985), respectively. A key paper on two-dimensional fields 
and its imphcations for the CMB is that of Bond & Efstathiou (1987). This technique 
was applied to calculate the number of spots on small angular scales in different models 
(Vittorio & Juszkiewicz 1987; Martinez- Gonzalez and Sanz 1989) and to study the Tenerife 
experiment (Gutierrez et al. 1994). A similar analysis applied to non-gaussian random 
fields have been performed by Coles & Barrow (1987). 

In this Paper, we consider the statistical properties of the CMB, assuming that the 
temperature fluctuations can be represented by a two-dimensional gaussian random fleld. 
The local description of maxima is presented in §2. We will restrict our analysis to the 
peaks of the field above a certain threshold. In particular we are interested in the following 
quantities: mean number of peaks over the whole celestial sphere A^(> u) (§3), gaussian 
curvature probability density function (p.d.f.) p(«;, > u) (§4) and eccentricity p.d.f. p(e, > v) 
(§5). All of them can be calculated analytically in terms of parameters that are related to 
the CMB radiation power spectrum. This power spectrum is characterized by the Doppler 
peaks and a cut-off at high £, which depend on the cosmological parameters (for a recent 
review see Hu 1996). An accurate determination of the radiation power spectrum requires 
detailed numerical calculations in perturbation theory (Sugiyama 1996). Flat and open cold 
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dark matter models with a Harrison-Zel'dovich initial power spectrum will be considered. 
Our analysis includes an angular resolution ranging from 5' to 20', of interest for future 
experiments, and also the effects of noise. Conclusions are presented in §6. 



The local description of maxima involves the second derivatives of the field along the two 
principal directions. As usual, the curvature radii are defined by: Ri — [— A'/(max)/2]~^/^ 
and i?2 = [— A2(maa;)/2]^^/^, where A is the temperature field normalized to the 
rms-fluctuations. Then with any maximum, we can associate two invariant quantities: the 
gaussian curvature k and the eccentricity e given by 



The number density of peaks of a two-dimensional homogeneous and isotropic gaussian 
random field has been studied by Longuet-Higgins (1957) and Bond and Efstathiou (1987). 
After a straightforward calculation, one can obtain the mean number of maxima (over 
the celestial sphere) N{K,e,i')dKdedi' with gaussian curvature, eccentricity and threshold 
between {k, k + dn), (e, e + de) and {u, v + du), respectively. N is given in terms of two 
spectral parameters 7 and 9^ that charaterize the intrinsic cosmological model plus the 
noise 



2. Local description of mcixima 




(1) 



,1/2^ 



7 




(2) 




C(0,a), al = -2C"{Q,a), 




(3) 



where a is the gaussian dispersion (cr = 0.425 x FWHM). The C^'s have two different 
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contributions: the intrinsic cosmological signal and the noise. We will consider flat and 
open CDM models {Vl = 1,0.3,0.1, baryon content Qf, = 0.05, Hubble constant h = 0.5) 
with adiabatic fluctuations and a Harrison-Zel'dovich primordial spectrum, kindly provided 
by N. Sugiyama. The C^'s have been normalized to the COBE 2-year maps (Cayon et al. 
1996. This normalization does not appreciably change with the 4-year data). We consider 
several noise amplitudes assuming that all scales contribute at the same level (white noise). 
The multipole coefficients of the noise C(jv) are therefore given by: 

47r 

C(iv) = ^?iv)(10')^^(2^^1)e-^(^+iV ■ 

The noise amphtude 74(jv)(10'), i.e. the noise after smoothing with a 10' FWHM 
gaussian window, is fixed at A(^j^){10') — (0,1,3) x 10~^ with a — 0.425 x 10', giving 
C(Ar) = (0, 1.9, 17) X 10^^"'' for the three noise levels used in our examples (a justification will 
be given below). Then, A(^n^ for other angular resolutions (5', 20') can be obtained using the 
same C^n) through the previous formula. 

Following standard observational procedures, we will filter signal plus noise with 
a gaussian with approximately the same width than the antenna FWHM. The angular 
correlation function C{a,a) is therefore given by 

C{a, a) = j- Y.{2i + 1) (Qe-^(^+^)'^' + C(^)) P,(cosa)e-^(^+^)-^ (5) 

In Table 1 we give the coherence length and the parameters 7 and ^* for Q = 0.1, 0.3, 1, 
different angular resolutions and yl(jv) (10') 's. 6** increases with beam size and decreases 
with 74(jv)(10'). When no noise is present the coherence angle has a range between 8'. 6 and 
35'. 9 for the values of the parameters considered. That range decreases as the noise level 
increases. 
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We will analyse 2D temperature fluctuations (signal plus noise) with angular resolution 
FWHM{') = 5, 10 and 20, which are of interest for the most sensitive bolometers and 
radiometers of future space experiments (COBRAS/SAMBA and MAP) and also for the 
VSA experiment as well as other interferometric experiments. The values of A(jv)(10') 
considered in this paper cover the range of sensitivities expected for the future experiments. 
In particular, the best expected sensitivity of COBRAS/SAMBA corresponds, in practice, 
to the case A(jv)(10') — 0. 



3. Number of peeiks 

The number of peaks above the threshold N{> u), can be calculated from the 
differential number N{i>)di' 



N{v) 



Nt\- 



TT 



-1/2/2 



1' 



(z/^ ^erfc(7z/s)^ + 



i/7 



(1- 



1' 



TT 



1/2- 



+ W 1 _ erfc(7z/st) e 



(6) 



s=[2{l-j%"\ (=(3-2f)-'^ (7) 

and Nt — (s^/^^*^) is the total number of peaks over the whole celestial sphere. 

In Figure 1, we show the cumulative number of peaks for different values of Vt and 
A(iv)(10'). Generically, the number of peaks increases if we decrease either the beam size or 
the n parameter, except for FWHM = 5' with A(jv)(10') = 3 x 10-^ and FWHM = 20' 
where the number is greater for Q, = 0.3. For a noiseless map (i.e. yl(7v)(10') = 0) and 
angular resolution of 5', the number of peaks above the threshold 1/ = 3 for 17 = 0.1 is 
approximately 3 times the value for Q = 1 (i.e. 4541 as compared to 1657 peaks for the 
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open and flat cases, respectively). However, when noise is present, the most favourable 
case is an angular resolution of 10', at which the noise decreases considerably while the 
signal is slightly affected. In fact, using a test and assuming Poissonian errors, the 
hypothesis that the flat and open models are derived from the same population is rejected 
at a confldence level ^ 99% except for an angular resolution of 5', A(^n){10') — 3 x 10~^ and 
fl = 0.3. Moreover, as previously indicated, the best confidence level is attained at the 10' 
angular resolution. Since we are considering very small angular scales, the cosmic variance 
will not affect our results from the practical point of view. 

On the other hand, we may ask whether gravitational lensing can change these 
results. An estimate of the coherence length including, or not including, lensing leads to 
{Of /9cY ~ 1 — a^, with a = {a{P)//3)p=o being the relative bending dispersion at zero lag. 
For standard CDM and low-il CDM models: a ^ 0.18, so the number of maxima (which 
at high thresholds is approximately proportional to the coherence length) is only slightly 
modified, i.e. ^ 3% (Martinez- Gonzalez, Sanz and Cayon 1996). 

In Table 2, we give the number of peaks above the thresholds u — 3, 3.5, 4 for different 
values of Q, angular resolutions and levels of A(jv)(10'). 

An equivalent quantity that can be used is the mean area of the peaks above a certain 
threshold (defined as the total area above that threshold divided by the corresponding 
number of peaks). The behaviour of this quantity can be easily obtained from the 
number of peaks and so it does not incorporate any new information that discriminates 
between the different models. As an example, for the case of FWHM=10', u — 3 and 
yl(Ar)(10') = 3 X 10~^, we find a mean area {arcmin^) of 42.4 for Q = 0.1 and 46.7 for Q=l, 
whereas for the 74(jv)(10') = these values increase to 132.0 and 266.1, respectively. 
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4. The distribution of gaussian curvature 

The distribution of peaks above the threshold v with inverse of the gaussian curvature 
L = between (L, L + dL), p{L, > u), can be obtained from the following p.d.f. 



= ( -y^'a^tL-^e"'^"'e-i*'^'erfc 



TT 



a = 27^^' (9) 



In Figure 2, we represent the p.d.f p{L, > u) for different angular resolutions 
(FWHM(') = 5, 10, 20) and threshold v = 3 . In all the cases, except when the beam size 
is very small (5') and the noise amplitude 74(jv)(10') high, the curves associated with flat 
and open models clearly differ. In those cases, using a KS-test, the null hypothesis that the 
flat and open models are derived from the same population is rejected at a confidence level 
^ 99%. The best case is obtained for an angular resolution of 10' when noise is present. 
Increasing the threshold slightly modifies the shape of the distribution: the height of the 
maximum increases and the curve is shifted to lower L, i.e., the peaks fall more rapidly for 
higher u . 

On the other hand, we can obtain the mean L for the different models. For the 
case of FWHM=10' and A^n){10') = 3 x 10-^ we find (in arcmin^) < L >= 36.0 and 
38.7 for fl = 0.1 and 1, respectively. If we consider the same cases for 74(jv)(10') = the 
corresponding mean L's are given by 106.9 {Q — 0.1) and 177.7 (Q = 1). Then, since the 
error in < L > due to cosmic variance is expected to be very small for the small angular 
scales considered, we can also use the mean values of L to distinguish between the fiat and 
open models. In Table 3 the mean L's are given for the models considered. 

In order to measure the gaussian curvature from a map obtained by an experiment, 
the required pixel size would need to be approximately one fifth of the typical curvature 
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radius of the maxima. The two curvature radii for each peak can be measured by a fit to a 
paraboloid centered on the maximum temperature. The pixel size should be a compromise 
between having an appropriate number of pixels to perform the fit and remaining in the 
vicinity of the maximum. In particular, if we want to test Q values as low as 0.1 for an 
angular resolution of 10' and ^(jv)(10') = the required size should be fa 2'. 



5. The distribution of eccentricities 

The distribution of peaks above the threshold v with eccentricity between (e, e + de), 
p(e, > I/), can be obtained from the following p.d.f. 



p(e,.) = ^^e-^^^-^^f^((i/.r)VV^(l-lerfc(i7VV.)) 

(5/737(1 - 7^)V + H\^vf{l - 7^))) , (10) 



H- (2-^^ C^H^^l^ (ID 

(3-272)e4 + 4(l-e2)' ^ ^(2-e2)2- ^^^) 

We note that there is an error in the expression given by Bond and Efstathiou(1987) 
for the conditional probability P(e|i/), where e is the ellipticity related to e by 
e = 2(e/(l + 2e)y^'^. We have studied p(e, > v) for different angular resolutions, noise 
amplitudes A(jv)(10'), thresholds and models. The main conclusion is that it would be 
difficult to distinguish between the cosmological models based on the comparison of 
eccentricities. As a typical example, in Figure 2, we represent the p.d.f p{e, > v) for the 
angular resolution of 10', threshold = 3 and no noise. The introduction of some level of 
noise clearly makes things worse. Hence we can generically say that the eccentricity p.d.f. 
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is not a good test for distinguishing between flat and open models. On the other hand, 
Gurzadyan & Kocharyan (1992, 1993) argue that mixing of photons in a universe with 
negative curvature will produce elongated shapes as compared to the flat case. All of our 
results lead to the opposite conclusion: the eccentricity p.d.f. for flat and low-fi universes 
show similar bell-shape (for thresholds u — 3,4:) with mean value < e >^ 0.7 and almost 
the same dispersion. Therefore, we conclude that the eccentricity is a bad discriminator of 
the fl parameter. 



6. Conclusions 

We have studied the distribution of peaks above a threshold u using the mean number 
and two local quantities: gaussian curvature and eccentricity. We have considered a whole 
sky coverage, with angular resolutions of 5, 10, 20 arcmin (antenna FWHM) and different 
levels of noise A(jv)(10') — (0,1,3) x 10~^, and we have calculated the distribution of 
these quantities for flat and open CDM models (with a Harrison-Zel'dovich primordial 
spectrum). Our main conclusions are that the number of peaks and the gaussian curvature 
are good discriminators of the geometry of the universe, whereas the eccentricity cannot be 
used to distinguish between different Q values. For thresholds u — 3,4, these curves are 
indistinguishable for flat and open models, and we disagree with Gurzadyan & Kocharyan 
(1992, 1993) who argue that mixing of photons in a space of negative curvature would tend 
to elongate the spots in the CMB. On the other hand, an angular resolution of 10' is the 
most appropiate to distinguish between low-Q and flat models when noise is present. 

We would like to thank N. Sugiyama for providing us the radiation power spectrum 
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support from the Spanish DGICYT, project PB92-0434-C02-02. RBB acknowledges a 
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Fig. 1. — Logarithm of the number of peaks above a threshold v versus the threshold, for 
different angular resolutions. The dashed and dotted lines correspond to open universes (Q 
=0.1 and 0.3, respectively) and the solid one to a flat universe. Each set of three different 
lines corresponds to values of A(iv)(10') = (3, 1, 0) x 10~^ (from top to bottom). 

Fig. 2. — We plot in three of the figures the distribution of L for peaks above a threshold 
V — ?> and several angular resolutions. The first, second and third set of 3 different hues 
(from left to right) corresponds to values of A{jq-^{\SS) = (3, 1,0) x 10~^, respectively. The 
p.d.f. p(e, > 1/ ) for a signal-dominated map (74(jv)(10') = 0), u =3 and FWHM=10' is shown 
in the bottom right figure. In the four plots the dashed, dotted and solid hues correspond 
to the cases Q =0.1,0.3 and 1, respectively. 



Table 1: Coherence angle and 7 and 6** parameters related to the peaks 
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Table 2: Number of peaks above the threshold v 
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Table 3: < L > (arcmir?) of the peaks above the threshold v = Z 

^(jV)(10') = A(Ar)(10') = 10-5 yl(Ar)(10') = 3 X 10-5 

FWHM 0=0.1 0=0.3 0=1 0=0.1 0=0.3 0=1 0=0.1 0=0.3 0=1 
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